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SUMMARY 


A capability for the nonlinear vibration analysis of beam and frame 
(structures suitable for use with NASTRAN Level 15.5 is described. The 
nonlinearity considered is due to the presence of axial loads induced by 
longitudinal end restraints and lateral displacements that are large com- 
pared to the beam height. This paper includes a brief discussion of the 
mathematical analysis and the geometrical stiffness matrix for a prismatic 
beam (BAR) element. Also included are a brief discussion of the equivalent 
linearization iterative process used to determine the nonlinear frequency, 
the required modifications to subroutines DBAR and XMPLBD of the NASTRAN code, 
and the appropriate DMAP ALTERS to determine the frequency. To demonstrate 
this nonlinear vibration capability, four example problems are presented. 
Comparisons with existing experimental and analytical results show that 
excellent accuracy is achieved with NASTRAN in all cases. 
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INTRODUCTION 


In practical beam vibration problems, transverse deflections may be 
greater than those assumed for linear theory. Considerable attention has 
been given, therefore, to the nonlinear flexural vibration of beams. Most 
studies have dealt with simple uniform beams with either hinged or fixed 
support conditions at both ends. An excellent literature survey is given 
by Eisley (ref. 1) through 1964. More recent surveys are given by Ray and 
Bert (ref. 2), and Pandalai (ref. 3). Nonlinear vibrational behavior of 
non-uniform beam and frame structures found in many engineering applications, 
however, has not received much attention in the literature because of 
analytical difficulties. The use of the finite-element method overcomes 
these difficulties and removes the uniform beam and limited support condition 
restrictions. 
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The purpose of this paper is to describe a nonlinear vibrational analysis 
capability for determining fundamental frequency of beam and frame structures 
suitable for use with NASTRAN Level 15.5 and to present results demonstrating 
this capability. The paper includes a discussion of the mathematical analysis 
and the derivation of the geometrical stiffness matrix that represents the 
induced axial force in the governing equation, the appropriate modifications 
to the NASTRAN code, and solutions of example problems. Procedures for non- 
linear vibration analysis with and without applied axial forces are available 
for NASTRAN Level 15.5 by means of DMAP ALTERS and modifications of the 
NASTRAN code given in the Appendices. 


SYMBOLS 
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E 

{f} 
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I 

[k] 

[k d ] 

[k 8 ] 
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[m] 


P 


xo 


p 

X 


[p] 

{u} 


area 

amplitude of vibration 
amplitude ratio, c/^ l/A 
modulus of elasticity 
element forces 
height of beam 

area moment of inertia of cross section 

stiffness matrix 

differential stiffness matrix 

geometrical stiffness matrix 

length 

mass matrix 

applied axial force on undeformed beam 
axial force due to deflection 

Euler load 

inplane force matrix 
nodal displacements 
time 
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x,y,z 

element coordinate system 

w 

lateral deflection 

P 

mass density 

{<}>} 

normal mode 

0) 

0 

fundamental linear frequency 

0) 

fundamental nonlinear frequency 

Subscripts : 


a,aa 

system 

e,ee 

element 

eq 

equivalent 


MATHEMATICAL ANALYSIS 
Formulation of Matrix Equation of Motion 


The strain energy in a deformed bar element of uniform cross section is 
given by: 



The first two terms of the strain energy expression are due to bending and 
applied axial force, respectively* The last term is the nonlinear contribution 
of the axial force P induced by large deflections. If the axial force P x is 
neglected from Eq. (1) , the strain energy is reduced to that of the 
linear theory. 


The kinetic energy is given by: 
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where the rotatory inertia term has been neglected. 


In the finite-element approach, transverse deflection of an element is 
expressed in terms of generalized displacements {u } and interpolation 
functions. The function chosen to represent the transverse deflection is 
given by: 

w = |^(1 - 3 + 2 x^) -(1 - 2 x + x^)x (3 x^ - 2 x^) 


(x - X 2 ) X j {u } 

J e 


where x = ^ 


This displacement function is the same as used for the NASTRAN bar element. 
The displacement vector describing bending of a bar element in the xz plane 
(see figure 1) is defined by: 


{u } = {u 0 U 0 - } 
e za ya zb yb 


Substituting Eq. (4) into Eqs. (1) and (3), gives: 


u = i (u } T ( [k ] + [k d ] + [k 8 ]) {u } 

Z e ee ee ee e 


T = \ {u } T [m ] {u } 
2 e ee e 


where [k ] , [k ], [k g ], and [m ] represent the element stiffness, 
ee ee ee ee 

differential stiffness, geometrical stiffness, and mass matrices, respectively 
and [k^ ] is a function of {u }. Substituting Eq, (4) into Eq. (2) yields 


P = {u } [p ] (u } 
x e ee J e 


where [p ] is the inplane force matrix. 
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Substituting Eqs. (7) and (8) into Lagrange's equations, that is. 




i - 1, 2, . . , n 


(10) 


where n is the number of elemental degrees of freedom, leads to the matrix 
equation of motion for the large amplitude free oscillations of a bar element 
which is given by 


[ml (u } + ([k ] + [k d 1 + [k g ]){u } = {f } 


ee 


ee 


ee 


ee 


(id 


,The mass, stiffness, and differential stiffness matrices of the bar element 
are 12 x 12 matrices relating the forces and moments acting at the ends of the 
bar (see ref. 4). The portion of the differential stiffness matrix, for 
example, that describes bending in the xz plane of figure 1, is given by 


u 

za 

0 

ya 

u zb 

6 ^ 
yb 

6 P 

P 

6 P 

P 

xo 

xo 

xo 

- xo 

51 

10 

51 

10 


(k° 


ee 


1 = 


2SL P P 

xo xo 

15 10 


l P 

xo 

30 


( 12 ) 


6 P 

P 

xo 

xo 

51 

10 


21 P. 

15 


The relations between [k d ] , [k 8 ], and [p ] can be found from Eqs. (1) 

and (2) and they are 6 “ e 6 


[k 8 1 = | [k d ] 

ee z P ee 


xo 


(13) 
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where 


P 

x 



ip ee > 


{u } 

e 


and 


(14) 


[p 


ee 



P 

xo 


tk d J 

ee 


(15) 


Solution Technique 

Q 

The goemetrical stiffness matrix [k ] in Eq. (11) is displacement 
dependent (see Eqs. (13) and (14)). Therefore, the frequency for 

nonlinear vibration also depends on the amplitude of vibration. This 
phenomenon is different from the linear case, in which the frequency is 
independent of amplitude. In the following the frequency associated with 
the linear vibration problem is referred to as the linear frequency, and the 
frequency associated with the nonlinear vibration problem is referred to as 
the nonlinear frequency. To determine the nonlinear frequency, an iterative 
procedure with an equivalent linearization technique is used and is illus- 
trated by the simplified flow chart shown in figure 2. The system matrices 
indicated in figure 2 are assembled from the element matrices by a standard 
finite-element procedure. The basic idea is to replace the displacement depen- 
dent geometrical stiffness matrix [k^ ] by an equivalent matrix [k^ ee ] e q 

using the mode shape of the linear vibration problem as the first 

approximation to the displacement. This reduces the nonlinear system 
equation of motion to a linearized equation which can be solved as a standard 
eigenvalue problem. The mode obtained by solving this eigenvalue problem 
may be used to recompute [k^^] for the next iteration in the nonlinear 
vibration iterative solution procedure. The solution procedure is illustrated 
as follows. The first step is to solve the linear vibration problem: 

0) 2 [m ] {<j>} = [k 1 {<}>} (16) 

o aa o aa o 


where cj is the fundamental frequency of the linear problem, {(f)} represents 
the corresponding mode shape normalized by its maximum component, and the sub- 
script aa represents the system matrices. Solving Eq. (16) provides the first 
approximate displacement in the form 



inhere c is the amplitude of vibration. The equivalent geometrical stiffness 
natrix now can be obtained using {u }]_ which leads to a linearized eigenvalue 
aquation of the form a 

i 

I 


w I“aa 3 { * } 1 “ ([k aa ] + [k aa^ + ^aaW^l 


(18) 


i 

fcrhere w is the fundamental nonlinear frequency associated with amplitude c, 
and is the corresponding mode shape. The iterative process can be 

repeated by using 


{u a>2 * c (19) 


as the second approximation, and similarly the i-th iteration approximate 
displacement is of the form 


{u }. = c 
a x 


Wi-i 


( 20 ) 


The iterative process can be continued until the nonlinear frequency converges 
to the desired accuracy or the mode shape {<}>}. satisfies some convergence 
criterion (e.g., the modified Euclidean norm 1 of ref. 5). 


MODIFICATIONS TO THE NASTRAN CODE 


To compute the geometrical stiffness matrix [k^ ], subroutine DBAR was 

modified to take advantage of the fact that the differential stiffness matrix 
and the geometrical stiffness matrix are related as shown in Eq. (13). 

Appendix A shows these changes in CDC UPDATE format. The core storage require- 
ment for DBAR was increased by 647 0 locations. 

o 

To avoid going through the modified section of code each time DBAR was 
| called, a new parameter, NLVIB, was added to the DMAP calling sequence for 
'module DSMG1. The contents of NLVIB are passed through blank common from 
DSMG1 to DBAR. The default value for NLVIB, set in block data subroutine 
XMPLBD, is zero (0). When NLVIB = 0, the new code in DBAR will not be 
executed. To set NLVIB = 1 and execute the new DBAR code, the following 
calling sequence for the DSMG1 module is used: 

DSMG1 CASECC, ,SIL, ,PHIG,CSTM,MPT,ECPT,GPCT, 

DIT/K GGG/V ,N,DSCOSET /C ,N, 1 $ 

The underlined parameter sets NLVIB to 1. 
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Once the changes shown in Appendix A were made to DBAR and XMPLBD , they j 
were compiled and replaced the old DBAR and XMPLBD in the NASTRAN object 
library. Link 1 and Link 13 were relinked and a new executable NASTRAN was 
created. Although this procedure was done on a CDC computer, similar proce- 
dures will produce similar results on both the IBM and UNIVAC computers. In 
order to use this nonlinear vibration capability in NASTRAN, extensive alters 
have to be applied to either Rigid Format 5 and Rigid Format 13, depending 
on how the capability is needed. A summary of the applicable analyses, their 
governing equations and their appropriate Rigid Formats is given as follows: 


ANALYSIS 


Normal Modes 
Buckling 


EQUATION 
w o 2 [m] = [k] 

[k] + A[k d ] = 0 


RIGID FORMAT 
3 


Normal Modes with 2 ^ 

Differential U) [m] = [k + k ] 13 

Stiffness 


Nonlinear Vibration 
Analysis 

Nonlinear Vibration 
Analysis with 
Differential 
Stiffness 


d) 2 [m] = [k + k 8 ] 5 with ALTERS 

co 2 [m] = [k + k d + k 8 ] 13 with ALTERS 


where X is an eigenvalue. 

The appropriate DMAP alter sequences for both Rigid Formats 5 and 13 are 
shown in Appendix B. The alters between the statements 

LABEL CONV $ 

REPT CONV, 1 $ 

will go through two iterations. If the user desires more iterations, the 
integer in the REPT statement must be increased. The only other input re- 
quired to use this capability is the addition of a PARAM card in the Bulk 
Data deck. The parameter AMP is used to specify the amplitude of vibration 
of the structure. 


266 


I 


EXAMPLES AND RESULTS 


The nonlinear vibration capability developed for use with NASTRAN has 
been demonstrated by solving two examples of a nonlinear vibration analysis 
and two examples of a nonlinear vibration analysis with differential stiff- 
ness. NASTRAN solutions are compared with previously published results. 


| Nonlinear Vibration Analysis 

I 

The first example is the vibration of a uniform beam with various end 
| support conditions. Evensen (ref. 6) obtained approximate amplitude- 
frequency relations for uniform beams with fixed-fixed, hinged-hinged, and 
. fixed-hinged boundary conditions using a perturbation method. Good agree- 
ment is obtained between the NASTRAN and perturbation solutions as shown in 
i Figure 3. For the hinged support case, the two amplitude— frequency curves 
I coincide. 

i 

The second example demonstrates the effect of the amplitude of vibration 
on a rectangular frame structure. The frame is 304.8 cm (10.0 ft.) wide, 
609.6 cm (20.0 ft.) long, and is made of 1.27 cm (0.5 in.) diameter steel rod. 
There are 10 equally spaced cells lengthwise and 4 equally spaced cells along 
I the width. All four edges of the frame are fixed. A plot of the undeformed 
frame is shown in figure 4. Only one-fourth of the frame is used in the 
analysis due to symmetry. The linear frequency and nonlinear frequencies for 
values of the amplitude c (see Eq. 17) up to 7.62 cm (3 in.) are as follows: 


Amplituda, c Frequency, Hertz 


cm 

(in. ) 


0 

0 

4.638 linear 

1.27 

0.5 

5.319 

2.54 

1.0 

6.830 

3.81 

1.5 

8.565 

5.08 

2.0 

10.295 

6.35 

2.5 

11.940 

7.62 

3.0 

13.478 


The results indicate that the amplitude has an important influence on 
the frequency. In this example, a 5.08 cm (2 in.) amplitude at the center * 
of the frame increases the fundamental nonlinear frequency to more than twice 
the linear one. 
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Nonlinear Vibration Analysis with Differential Stiffness 


The third example is a hinged rectangular beam subjected to an applied 
axial tensile force of 105.4 N (23.7 lbf ) . The same problem was solved in 
ref. 2 using three different approximate analytical procedures, and results , 
from these procedures, as well as experiment results, are given. One procedure 
is based on an assumption for the spatial dependence of the displacement 
function, one is based on an assumption for the temporal dependence of the 
displacement function, and the third procedure is the Galerkin procedure. The 
beam has the following properties: 


50.8 cm (20.0 in.) 

1.27 cm ( 0.50 in.) 

0.081 cm ( 0.032 in.) 

Titanium Alloy r 

100.6 GPa ~ (14.6 x 10 psi) 

5.15 Mg/m ( 0.186 lb. /in. 3) 


Length 

Width 

Height 

Material 

Elastic Modulus 

Specific Mass 


One-half of the beam modeled by six BAR elements was used for the analysis 
Three analytical fundamental frequencies and an experimentally measured one fro 
ref. 2 and the NASTRAN solution are shown in figure 5, Comparing the results 
demonstrates that the NASTRAN results provide the closest comparison with the 
experiment . 

The fourth example is a beam-column subjected to an applied compressive 
force with various support conditions. Based on linear theory, Lurie (ref. 7) 
has shown that the relation between the square of the frequency and the axial 
load is linear for a beam that has identically shaped vibration and buckling 
modes. He also showed that the condition of zero fundamental frequency cor- 
responds to buckling. The linear vibration-stability problem studied by Lurie 
is actually the limiting case of a more general phenomenon of large amplitude 
vibrations under the influence of axial loads. Burgreen (ref. 8) obtained 
an exact solution in terms of elliptic functions for a uniform beam hinged at 
both ends and also verified his results experimentally. Srinivasan (ref. 9) 
used Galerkin' s method to study beam-columns with both ends hinged. Table 1 
shows good agreement between the NASTRAN solutions and the results given by 
Burgreen and Srinivasan, and NASTRAN gives better predictions than the one-term 
Galerkin method. The load-frequency curves for different amplitude ratios of 
vibration, d, for beams with various support conditions are given in figure 6, 
where d is the ratio of amplitude to the radius of gyration of the beam. No 
comparison is made for the cases of fixed-hinged and fixed-fixed because there 
appears to be no solution available in the literature. It has been found from 
this example that (1) the effect of amplitude is more pronounced for a less 
stiff structure, and (2) nonlinear theory shows that the frequency of a column 
at the Euler buckling load is not zero for finite amplitudes of vibration. 
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CONCLUDING REMARKS 


Nonlinear vibration capability for beam and frame structures has been 
developed for use with NASTRAN Level 15.5 by means of DMAP ALTERS and modi- 
fications to the NASTRAN code. A geometrical stiffness matrix for a bar 
element has been developed for NASTRAN by modifying subroutine DBAR. An 
equivalent linearization technique and iterative process used to determine 
nonlinear frequencies are implemented into NASTRAN by the DMAP ALTERS. The 
versatility of the finite-element method enables the analyst to determine 
nonlinear frequencies of vibration for non-uniform beam and frame structures. 
Comparison with previously published results show that excellent accuracy 
| is achieved with NASTRAN. 

i 
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APPENDIX A 


MODIFICATION OF CODE 


* I NSERT « DBAR • 1 1 1 

DOUBLE PRECISION PXT « COEFF » HOLD ( 6 ) « TEMP ( 36 ) < PXP « STOP ( 36 ) 
*INSERT»OBAR. 185 
C 

C INSERT NEW NL V 1 6 PARAMETER INTO oLANK COMMON 
C 

COMMON I COM, NLV IB 

* I NSERT .DBAR.638 

IF (NLV IB.NE.O# AND.FX.EQ.O.ODO ) FX=1 .ODO 
*INSFRT , DBAR. 646 

* I NSERT . DBAR . 766 
IF (NLV IB.EQ.O ) GO TO 621 

DIVIDE FX OUT OF KDGG MATRIX 
C 

DO 900 KK = 1 i 144 
KD (KK ) = KD ( KK ) /FX 
900 CONTINUE 
C 

C COMPUTE 

C T E*A 

C PXP = U * * (KDGG) * U 

C 2*L 

C 

C WHERE U is a 6X1 VECTOR 
C 

C KDGG IS A 12X12 MATRIX DIVIDE INTO FOUR 6X6 MATRICES 

C 

IF ( IPVT.NE. 1 ) GO TO 62 1 
COEFF=DA*E/(2.0D0*L ) 

PXT = O.ODO 
DO 920 IP1 = 1 .4 

JPX = 1 

JCNT = 0 

I F ( I P 1 . GE . 3 ) GO TO 930 
I LO= 1 
I H I =72 
GO TO 940 
930 I LO=73 
I H I = 1 44 

940 DO 960 I = I LO . I H I » 1 2 
JLO= I 
JH I = JLO+5 
DO 950 K = JL0.JH1 
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APPENDIX A 


I 

1 JCNT = JCNT+1 

I F ( I P 1 • EQ • 1 .0R.IP1.EQ.3) TEMP (JCNT ) =KD(K ) 
j IF< IP1 .EQ.2.0R. IP1 .EQ.4 ) TEMP < JCNT ) =KD ( K+6 ) 

1950 CONTINUE 
960 CONTINUE 
I WLEFT = 37 

I F ( I P 1 • GE • 3 ) I WLEFT = 73 
I WRGHT = 37 

IF< IP1 .EQ.2.0R. IP1 .EQ.4 ) IWRGHT=73 

CALL GMMATD ( KE ( I WLEFT >.6,6.1 . TEMP ( 1 ) ,6,6,0.3TOR( 1 ) ) 

CALL GMMATD < S TOR ( 1 ) , 6,6,0 » KE ( I WRGHT ) ,6,6,0, TEMP ( 1 ) ) 

I 912 DO 915 I P2 = 1 ,36 

TEMP ( I P2 ) = TEMP < IP2 )*COEFF 

|915 CONTINUE 

i IF( IP1 .EQ.2.0R. IP1 .EQ.4 ) jPx = 2 

, I F ( I P 1 • LE .2 ) CALL GMMATD< UA < 1 ) ,6, 1 , 1 , TEMP ( 1 ) ,6, 6.0, HOLD ( 1 > ) 

> I F ( I P 1 • GE . 3 ) CALL GMMATD(UB ( 1 > , 6, 1 , 1 ,TEMP ( 1 ) ,6,6, O.HOLD ( 1 > > 

I F ( JPX .EQ.l) CALL GMMATD (HOLD( 1 ) , 1 ,6. O.UA ( 1 ) , 6, 1 , O.PXP ) 

I F ( JPX • EQ . 2 ) CALL GMMATD(HOL.D( 1 ), 1 ,6,0, UB ( 1 > ,6, 1 , O.PXP) 

0 

C STORE SUM INTO PXT 

f' 

PXT = PXT+PXP 
920 CONTINUE 
C. 

C PXT = • 5 * PXP * KDGG 
C 

PXT = PXT * . 5D0 
621 CONTINUE 
* 1 NSERT , DBAR .790 

IF (NLV I3.EQ.0 ) GO TO 653 
C 

C CALCULATE KDGGG AND STORE 
C 

DO 652 IPX=1 ,36 

JPX= IPX+1 08 

KEP ( JPX ) =PXT*KEP ( JPX ) 

| 652 .CONTINUE 
653 CONTINUE 
-X-COMPIL^.DBAR 
* DELETE , XMPLBD . 78 

1, 10, 4HDSMG , 4H 1 , 1,10, 1, 1, -1, 1, O 

*DELETE . WT 1 33-L 14.39 

DIMENSION MPL01 ( 48), MPL02(178), MPL03(191), MPL04 ( 1 79 ) 
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APPENDIX A 


♦DELETE . WT 1 55-L 16.8 

COMMON/ XGP I 2/LMPL » MPLPNT . MPL (1621 ) 

♦ DELETE * WT 1 33-L1 4.42 . WT1 33-L 14.45 

3 , (MPL< 227 ) , MPL03 ( 1 ) ) . (MPL< 4 18). MPL04 ( 1 ) ) 

5 » ( MPl_ ( 597 ) . MPL05 ( 1 ) ) . (MPL( 723 ) . MPL06 ( 1 ) ) 

7 . ( MPL < 858 ) . MPL07 ( 1 ) ) , ( MPL ( 1 067 ) . MPL08 ( 1 ) ) 

9 « ( MPL ( 1 1 72 ) . MPL09 ( 1 ) ) » ( MPL ( 1 35 1 ) » MPL 1 0 ( 1 ) ) 

♦DELETE . WT 1 55-L 16.9 

1 . ( MPL ( 1 538 ) ♦ MPL 11(1)) 

♦ COMP 1 LE « XMPL8D 
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APPENDIX B 


DMAP ALTERS 
For Nonlinear Analysis 


I D NLVBF.BEAM frame 
APP DISP 
SOL 5*0 
TIMF 5 

jj NONLINEAR FREE VIBRATIONS OF BEAM AND FRAME STRUCTURES 

ALTER 19* 23 

GP3 GE0 iv, 3»EQEXIN*GE0M2/*GPTT/C»N*123/V*N*N0GRAV/C*N»123 $ 

ICHKPNT GPTT $ 

ALTER 31*31 

[ALTER 54 

EQUIV MGG.MNN/MPCF 1 $ 

[CHKPNT MNN « 

ALTER 66*76 

|MCE2 USET,GM»KGG*MGG»*/KNN*MNN,, $ 

'CHKPNT KNN * MNN $ 

LABEL LBL? S 

EQUIV KNN*KFF/SINGLE /MNN ♦ MFF/S I NGLE * 

CHKPNT KFF ♦ MFF * 

COND LBL3*SINGLE % 

| SC E 1 USET*KNN*MNN* */KFF*KFS*KSS*MFF» ♦ S 

ICHKPNT KFS*KSS»KFF*MFF $ 

.LABEL LBL3 * 

EQUIV KFF ♦ K A A/OM I T /MFF ♦ MAA/OM I T $ 

CHKPNT KAA*MAA $ 

ALTER 79 

SMP2 USET*GO*MFF/MAA $ 

CHKPNT MAA S 

i ALTER 61*129 

ALTER 134,141 

jSETVAL //V,N*bREAK/C,N* 1 / V ♦ N * L I Nk/C , N * - 1 % 

SAVE BREAK »L INK $ 

LABEL CONV $ 

| EQU I V KAA*KDAA/BREAK $ 

[EQUIV MAA*MDAA/8REAK $ 

! READ KAA*MAA* , * EED , USET * C ASECC/L AM A * PH I A * M I ,OE I GS/C ♦ N , MODES/ V * 

N*NEIG/C»N*2 S 
| SAVE NEIG % 

CHKPNT LAMA,PHIA*MI ,0EIGS S 

OF P LAMA *OEIGS, * ♦ , //V * N , C ARDNO $ 

SAVE CARDNO $ 

COND FINIS* NEIG $ 

SDR 1 USET* *PHIA, , ,GO*GM* *KFS* */PHIG* *BQG/C,N, 1/C»N»REIG $ 

CHKPNT PHIG*BQG % 

EQUIV PHIG, PHI AMP/BREAK $ 

ADD PHIG*/PHIAMP/V*Y, AMP $ 
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APPENDIX B 


CHKPNT 

EQUIV 

DSM61 

CHKPNT 

ADO 

CHKPNT 

EQUIV 

CHKPNT 

COND 

MCE? 

CHKPNT 

LABEL 

EQUIV 

CHKPNT 

COND 

SCE 1 

CHKPNT 

LABEL 

EQUIV 

CHKPNT 

COND 

SMP1 

CHKPNT 

SMP? 

CHKPNT 

LABEL 

EQUIV 

EQUIV 

REPT 

ADD 

ADO 

alter 

endalter 

CENO 


PH I AMP S 

PHIAMP,PHIG/LINK $ 

C ASECC « ♦ S IL » *PHIG»CSTM*MPT ♦ ECPT ♦ 6PC T ♦DIT/KGGG/V»N* 
DSCOSET/C»N» 1 * 

KGGG $ 

KGGGfKGG/KDGGG S 
KDG6G $ 

KDGGG»KDNN/MPCF2 /MGG .MDNN/MPCF2 % 

KDNM ♦ MDNN $ 

LBL2D * MPCF2 * 

USE T « GM » K DGGG ♦ MGG * */KDNN» MDNN ♦ * t. 

KDNN ♦ MDNN S 
LBL2D % 

kdnn«kdff/single /mdnn,mdff/single 

KDFF ♦ MUFF $ 

LBL3D»SINGLE « 

USET ♦ KDNN ♦ MDNN ♦ ♦/KDFF,KDFS» * MDFF * , * 

KDFF »KDFS»MDFF $ 

LBL3D i 

KDFF»KDAA/OMIT /MDFF ♦ MD A A/OM I T $ 

KDAA»MOAA $ 

LBLSDiOMIT * 

USET » KDFF * ♦ ♦ /GDO»KDAA»KDOO*LDOO*UDOO* ♦ t » ♦ % 
GDOtKDAA $ 

USET ♦GOO»MDFF/MDAA $ 

MD A A % 

lblsd % 

KDAA *KAA/LINK % 

MDAA.MAA/LINK % 

CON V * 1 $ 

MDAA »KDAA/KMAA $ 

PHI A MP » /PH I M % 

1S7* lSd 
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Nonlinear Vibration Analysis With Differential Stiffness 


ID NLVDS.BEAM FRAME 

»PP DISP 

>01 13.0 

'IMF. 5 

> NONLINEAR VIBRATION MODE WITH DIFFERENTIAL STIFFNESS FOR BEAM 

» AND FRAmE STRUCTURES 

kLTER SO 

, Q U I V MGG ♦ MNN/MPCF 1 <6 

:hkpnt mnn * 

kL TER 62. 72 

1CE2 USET, GM.KGG.MGG. ./KNN.MNN, , % 

CHKPNT KNN.MNN $ 

LABEL LBL2 S 

EQUIV KNN.KFF/SINGLE /MNN . MFF/S I NGLE $ 

kHKPNT KFF.MFF % 

tOND LBL3 ♦ SINGLE f. 

5CE1 USET.KNN.MNN. ,/KFF.KFS.KSS.MFF,, % 

;hkpnt kfs.kss.kff.mff $ 

1.ABEL L BL 3 S 

•[QUIV KFF.KAA/OMIT /MFF ♦ MA A/OM I T $ 

tHKPNT KAA.MAA $ 

ALTER 7b 

SMP2 USET.GO.MFF/MAA $ 

DHKPNT MAA $ 

ALTER 106.12b 

ALTER 130.130 

SETVAL //V.N.BREAK/C.N. l/V.N, LINK/C, N,-l $ 

SAVE BREAK .LINK $ 

Label conv % 

EQUIV KAA.KSAA/BREAK $ 

EQUIV MAA.MSAA/8REAK $ 

READ KAA.MAA. , , EED , USET ♦ C ASECC /L AMA , PH I A ♦ M I . OE I GS/C , N . MODES/ V ♦ 

N.NFIGV/C.N.2 S 
ALTER 136.137 

SDR 1 USET, ,PHIA, , .GO.GM. ,KFS» ♦ /PH I G. ♦ BQG/C . N ♦ 1 /C , N , BKL 1 $ 

CHKPNT PHIG.BQG $ 

EQUIV PHIG, PHI AMP/BREAK $ 

ADD PHIG. /PHIAMP/V.Y. AMP S 

CHKPNT PH I AMP $ 

EQUIV PHIAMP.PHIG/LINK $ 

DSMG1 C ASECC, ,SIL, . PH I G , CSTM . MPT , ECPT . GPCT , D I T/KGGG/ V , N . 

DSCOSET/C.N, 1 % 

CHKPNT KGGG $ 

ADD5 KGG ♦ KDGG , KGGG » ,/KSGG $ 

CHKPNT KSGG S 
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APPENDIX B 


EQUIV 

KSGG* KSNN/MPCF2 / MGG ♦ MSNN/MPCF2 $ 

CHKPNT 

KSNN * msnn $ 

COND 

LBL2S*MPCF2 $ 

MCE2 

USET *GM*KSGG*MGG* »/KSNN*MSNN* ♦ 4 

CHKPNT 

K SNN » MSNN $ 

LABEL 

LBL2S $ 

EQUIV 

KSNN.KSFF/SINGLE / MSNN*MSFF/SINGLE 

CHKPNT 

K SFF * MSFF $ 

COND 

LBL3S»SINGLE 4 

SCE I 

USET* KSNN* MSNN* * /KSFF *KSFS ♦ ♦ MSFF ♦ ♦ 4 

CHKPNT 

KSFF * KSFS ♦ MSFF $ 

LABEL 

LBL3S 4 

EQUIV 

KSFF*KSAA/OMIT / MSFF * MSAA/OM I T 4 

CHKPNT 

KSAA*MSAA $ 

COND 

LBLBS*OMIT 4 

SMP1 

USE T, KSFF* » */GSO*KSAA,KSOO*LSOO*USOO 

CHKPNT 

GSO » KS A A % 

SMP2 

USET*GSO»MSFF/MSAA i 

CHKPNT 

MS A A $ 

LABEL 

L8LBS 4 

EQUIV 

KSAA *KAA/L IimK $ 

EQUIV 

MSA A *MA A/L I NK $ 

PEPT 

COM V » 1 S 

ADD 

MSAA*KSAA/KMAA 4 

ADD 

PHI AMP * /PH I M l 

enqalter 


CEND 
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TABLE 1. FREQUENCY RATIO OF A HINGED BEAM 


Axial Load, 
•^xo 

P E 

Method 

Amplitude Ratio, 



d = 7T7 a 


0.0 

1.0 

2.0 

3.0 


NAS TRAN 

1.0000 

— 

1.0889 

1.3183 

1.6260 

0.0 

Ref. 8 

1.0000 

1.0892 

1.3178 

1.6257 


Ref. 9 

1.0000 

1.0897 

1.3229 

1.6394 


NAS TRAN 

.8944 

.9928 

1.2401 

1.5631 

0.2 

Ref. 8 

.8944 

.9930 

1.2389 

1.5618 


Ref. 9 

.8944 

.9937 

1.2450 

1.5772 


NAS TRAN 

. 7746 

.8864 

1.1566 

. 

1.4976 

0.4 

Ref. 8 

. 7746 

.8864 

1.1543 

1.4949 


Ref . 9 

.7746 

.8874 

1.1619 

1.5125 


NASTRAN 

.6325 

.7653 

1.0666 

1.4291 

0.6 

Ref . 8 

.6325 

.7649 

1.0627 

1.4246 


Ref. 9 

.6325 

.7665 

1.0724 

1.4448 


NASTRAN 

.4472 

.6210 

.9682 

1.3570 

O 

00 

Ref. 8 

.4472 

.6194 

.9617 

1.3502 


Ref. 9 

.4472 

.6225 

.9747 

1.3739 


NASTRAN 

.0026 

.4309 

.8586 

1.2810 

1.0 

Ref . 8 

.0000 

.4236 

.8472 

1.2708 


Ref. 9 

.0000 

.4330 

.8660 

1.2990 


NASTRAN 



.7329 

1.2000 

1.2 

Ref. 8 



.7105 

1.1851 


Ref. 9 



.7416 

1.2196 
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Figure 3. Amplitude-frequency curves for various support conditions. 
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